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We study the nature of the confinement phase transition in d = 3 + 1 dimensions in various 
non-abehan gauge theories with the approach put forward in [l[. We compute an order-parameter 
potential associated with the Polyakov loop from the knowledge of full 2-point correlation functions. 
For SU{N) with A'^ = 3, . . . , 12 and Sp{2) we find a first-order phase transition in agreement with 
general expectations. Moreover our study suggests that the phase transition in E{7) Yang-Mills 
theory also is of first order. We find that it is weaker than for SU{N). We show that this can 
be understood in terms of the eigenvalue distribution of the order parameter potential close to the 
phase transition. 
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I. INTRODUCTION 

Understanding the confinement of gluons is a chal- 
lenging problem, because characteristic and significant 
quantities providing both analytical and numerical ac- 
cess are difficult to identify straightforwardly. In recent 
years, the study of the confinement mechanism based on 
the IR behavior of gauge-dependent correlation functions 
has turned out to be both fruitful and inspiring. Most 
prominently, this has lead to two different confinement 
scenarios: the Gribov-Zwanziger [2], Q scenario on one 
hand and the Kugo-Ojima scenario on the other hand 
[J|. This scenarios are mutually related and connect sig- 
natures of color confinement with the low-momentum be- 
havior of gluon and ghost correlation functions. They 
have been intensely investigated by a variety of nonper- 
turbative field theoretical tools such as lattice gauge the- 
ory as well as functional methods j^-Q- 

Even if color confinement is eventually properly ac- 
counted for by a corresponding IR behavior of correla- 
tion functions, indicating, e.g, positivity violation and 
the absence of colored states in the physical state space, 
the remaining pressing open question is the relation of 
color confinement to quark confinement. Typical quark- 
confinement criteria such as those based on the Wilson- 
loop or Polyakov-loop expectation value [l^l in quenched 
QCD have long remained inaccessible from the pure 
knowledge of low-order correlation functions of the gauge 
sector. For instance, a direct computation of the heavy- 
quark potential requires additional knowledge, e.g., of the 
quark-gluon vertex |ll| or of the static gluon correlator 
in the Coulomb gauge [l^l ■ 

In m a more direct relation between color confinement 
and quark confinement has been established. There, a 
thermodynamic order parameter, namely the Polyakov 
loop, has been computed from the Yang-Mills propa- 
gators, thereby establishing for the first time a direct 
link between the IR behavior of correlation functions 
and quark confinement. Aside from the above relation 



of color confinement and quark confinement, the rela- 
tion between quark confinement and chiral symmetry 
breaking has more recently been studied with the aid 
of so-called dual observables [Ij, [ij] . These formal and 
computational advances better our understanding of the 
mechanisms underlying confinement and chiral symme- 
try breaking and pave the way towards a first principle 
access to QCD with functional methods. 

In particular, the approach put forward in [l| gives us 
access not only to an order parameter for confinement, 
the Polyakov loop, but also to its full effective potential. 
The latter is a crucial input in Polyakov loop extended 
effective models such as the PNJL and the PQM mod- 
els [l^. In these models, the Polyakov loop potential is 
an external input the parameter of which are fixed to re- 
produce pure Yang-Mills lattice results. Evidently, these 
physics constraints do not completely fix the potential, 
leave aside its extension to full QCD. Different potentials 
have been studied, for a comparison see e.g. [la|, and the 
physics at finite chemical potential is in fact very sen- 
sitive to parameter changes. In this regard the present 
approach provides the opportunity for a qualitative im- 
provement of the above models since it allows to fix the 
Yan g-M ills potentials completely from first principles, see 
also 171. This has motivated further studies such as the 



delevopment of a framework for the study of QCD with 
two colors [I3 . Moreover the present approach underlies 
the fully dynamical continuum study of two-fiavor QCD 
at finite temperature and quark chemical potential put 
forward in |14| . 

The present work is devoted to a more detailed analysis 
of the interrelation between color confinement and quark 
confinement as deduced in [l| . In addition to providing a 
more comprehensive technical insight into the underlying 
ideas, we illustrate this interrelation and the generality 
of the approach by applying it to a variety of non-abelian 
gauge theories near the deconfinement phase transition. 
In fact, while the phase transition is of second order in 
SU{2) Yang-Mills theory, it is well-known from lattice 



simulations that a first-order phase transition occurs in 
SU{N) gauge theories with TV > 3. This brings up the 
question on how the nature of the phase transition is re- 
lated to the properties of the underlying gauge group. It 
has been conjectured in [19[ based on the order-disorder 
nature of the deconfinement phase transition with respect 
to center symmetry that the phase transition of SU{2) 
Yang-Mills theory should fall into the Ising universality 
class. This observation based on the symmetry proper- 
ties of the center of the gauge group does not necessarily 
extend to the other gauge groups with a center symmetry 
agreeing with a 2nd order universality class. In J20l.l2l| it 
has been conjectured that the dynamics near the critical 
temperature is sensitive to the mismatch of the num- 
ber of dynamical degrees of freedom in the confined and 
deconfincd phase. In accordance with this conjecture a 
first-order phase transition has been found for Sp(2) even 
though the center of the group is Z{2) [20, [21|. In the 
present paper we study 5*^(2) and E{7) gauge theory (the 
center of both groups is Z{2)) and compare the results 
to our findings for SU{N) gauge theories also in order to 
shed more light on the conjecture put forward in [20, [2l[ . 
Our study of the deconfinement phase transition is 
based on an order parameter related to the Polyakov loop 
variable, 
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more precisely on the expectation value of (Aq) in 
Polyakov gauge, see also ^^. Then, (Aq) is sensitive to 
topological defects related to confinement [22|, and also 
serves as a deconfinement order parameter. 

The effective potential of (Aq) is accessible from the 
knowledge of gauge correlation functions by means of the 
functional renormalization group (RG), [l[. As an addi- 
tional characteristic ingredient for a quantitative under- 
standing of the phase transition we introduce and iden- 
tify eigenvalue distributions of the order parameter which 
exhibit characteristic traces and facilitate a quantitative 
understanding of the behavior of the corresponding effec- 
tive potential. 

The paper is organized as follows; In Sect. Ulland lllll 
wc discuss general aspects of functional flows for a study 
of non-abelian gauge theories. In Sect. II VI we discuss how 
background-field RG fiows can be constructed from RG 
flows in Landau-gauge Yang-Mills theories. We discuss 
a sufficient confinement criterion in Sect. |V] before we 
present our study of the nature of the phase transition in 
SU{N), Sp{2) and E{7) gauge theory in Sect.|Vll 



II. FUNCTIONAL FLOWS AND 
OPTIMIZATION 

For our study of the phase transition of non-abelian 
jauge theories we employ the functional RG for the ef- 



fective action F^ [23|. This allows us to interpolate be- 
tween the initial UV action related to the classical ac- 
tion Tk=A — S and the full quantum effective action 
F = Ffe=o, being the IPI generating functional. The 
infrared (IR) regulator scale k separates the fluctuations 
with momenta p^ > fc^ which are already included in F^ , 
from those with smaller momenta which still have to be 
integrated out. The full RG trajectory is given by the 
solution to the Wetterich equation [t ~ ln(fc/A)), 
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where F^ ' denotes the second functional derivative with 
respect to the dynamical field $, collectively summariz- 
ing gluon and ghost fields in the present context. The 
super trace STr sums over momenta, internal indices and 
species of fields and includes a negative sign for the ghost 
fields. The regulator function Rk specifies the details of 
the Wilsonian momentum-shell integrations. See p4J - [27| 
for reviews on gauge theories. 

In the present context, we are interested in the ef- 
fective potential for an order-parameter quantity which 
is related to the local part of the full effective action 
F = Tk=o. The latter can formally be obtained from the 
integrated fiow, 
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Eq. ^ is an equation for the full quantum effective action 
F and has a priori no A dependence. A partial integration 
leads to 



F= i STr In F (2) 
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Note that the first term on the right-hand side does not 
depend on the regulator function. The initial conditions 
of the flow at A including possible subtractions are com- 
prised in the second line of Eq. Q . The second term in 
the first line of (|4]) reads 
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For general regulators, ATk is only finite upon subtrac- 
tions contained in the second line of Eq. Q, as dtTj^ 
does not vanish for large momenta. Whereas the repre- 
sentation (j3]) thus is of limited practical use in the general 
case, it is ideally suited for the determination of an order- 
parameter potential which is UV finite from the begin- 
ning, as is the Weiss potential. In fact, our quantitative 
results for the order-parameter potential are dominated 
already by the first term of Eq. (|4]) . For the remainder of 
this section, we will, however, be concerned with an op- 
timized strategy to evaluate the contributions from the 
AFfc term. 



The evaluation of this 'RG-improvemcnt term' ~ dtVk 
requires two nontrivial ingredients: full gluon and ghost 
propagators in the presence of an IR regulator, Gk = 
(Fj. + Rk)^^ and the flow of the inverse propagator 

(2) 

QiKk-, ■ This information has been made available in 
[9|, |28| where the full momentum dependence of Landau- 
gauge propagators has been computed within optimized 
RG flows. For earlier RG calculations, acquiring partial 
knowledge about Yang-Mills propagators and providing 
evidence for the Kugo- Oj ima / Gribov- Z wanziger confine- 
ment scenarios, see [sl. M. l29l. |3G| . 

Optimization of the RG flow has not only the advan- 
tage of a more stable and faster convergent numerical 
scheme; in a given truncation, it can actually be posed 
as a stability/convergence problem [31|. It also provides 
for a link to using propagators obtained from lattice sim- 
ulations, see below. Full functional optimization can be 
reformulated as the quest for a minimal flow trajectory 
for general functional flows [23]. In other words, for a 
given gap l/fc^ff of the propagator which constitutes the 
inverse of the physical infrared cut-off, the integrated op- 
timized flow is already as close as possible to the full the- 
ory, as the remaining flow trajectory is minimal. This 
results in the following pro pagator for the corresponding 
optimal regulator -Ropt |25| : 
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with fc^^ff = T'o'^'ik''). The propagator (F^'^ + i?opt)"^ 

(2) 

is already the full propagator for all eigenvalues of Fq ' 
belonging to specfFg } > fc^fj, and is identical to the gap 
for the remaining eigenvalues, specjFQ '} < k^ff. The 
choice ^ requires non-trivial field redefinitions. For the 

1/2 " 

gauge field we have A = Zj( A with dtA = and 
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where the right-hand side is evaluated at $ = 0. Eq. ([7]) 
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ensures 9tF^ ^(^off ~ ^q ) = 0, and hence 

ir['^-T^,'^)9ik!^-4'^)^o. 



(8) 



The conditions ^ and ([7]) allow us to provide the optimal 
regulator in an explicit form 25|: 
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With the choice ©, the flow of Green's functions can 
be computed within an iteration of the integrated flow 

(2) 

starting from an initial value for Fq . 



Let us elucidate aspects of the optimized flow in the 
context of the integrated flow (|3]) . For the optimized flow 
the relation Q follows from a direct integration of the 
flow: the flrst term in Eq. Q relates to integrating the 

the second term 



dfk^g contributions of the related flow 



(2) 



in Eq. ^ is the t integral of the contributions ^ dfT^ 
Details of the numerical computation of Eq. Q can be 
found in App. [X] 

From a general perspective, the effective action ^ to- 
gether with the optimized regulator can be considered as 
a DSE within a consistent BPHZ-typc non-pcrturbative 
renormalization [25|, where the A-dependcnt terms pro- 
vide the classical action and the subtraction terms. The 
computational benefit in comparison to standard DSE 
equations is the explicit finiteness of Eq. 1^ in any trun- 
cation without the need of further additive or multiplica- 
tive renormalizations. The second term on the right-hand 
side of the first line constitutes an RG improvement term. 



III. BACKGROUND FIELD FLOWS 

In order to arrive at the effective potential for an order- 
parameter field, we parameterize the fluctuations with 
respect to a background field which is related to the or- 
der parameter. In Yang-Mills theories, this decomposi- 
tion into fluctuating modes and the background fleld can 
be organized such that the resulting background-field ac- 
tion preserves a residual gauge symmetry, e. g. |32| . This 
approach using the background-field gauge can be under- 
stood as a simple extension of Yang-Mills theories within 
general covariant gauges. The gauge condition 9^^^ ~ 



is generalized to 



D^{A){A-A)^^0, 



(10) 



for an unspecified background field A. Equation (|T0)) 
implemented on configuration space in a strict sense de- 
fines Landau-DeWitt gauge. A less strict Gaufiian aver- 
age over the gauge condition D^{A){A — A)fj = C with a 
probability distribution ^ exp(— 1/^ J trC^) leads to the 
background-field equivalent of a general covariant gauge. 
Such a formulation has the benefit of an auxiliary gauge 
symmetry for the effective action under a transforma- 
tion of both, the full gauge field A — ;> A + Du! and the 
background A ^ A + Dlo. In this manner, the gauge 
condition (jlOp is unchanged since a = A — A transforms 
as a tensor, a -^ [a,uj]. 

With the gauge fixing ^TU\\ . the effective action now 
depends also on the auxiliary field, F = F[$,A] with 
$ = {a,C,C). We emphasize that the background-field 
gauge transformation is an auxiliary symmetry. The ef- 
fective action F = F[$,A] still carries non-trivial sym- 
metry constraints, namely the Slavnov- Taylor identities 
(STI). These follow from a gauge or BRST transforma- 
tion of the field $ at fixed background. Indeed, the 
underlying STI are that of a standard covariant gauge. 
Even though the gauge invariance is an auxiliary symme- 



try, it facilitates the construction of a (physically) gauge- 
invariant effective action T[A] = r[0,A = A]. The flow 
equation for r[$, A] in such a setting reads 



d,T,[^,A]^lTr-^-^ 
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The action rfc[$,A] is still gauge invariant under back- 
ground gauge transformations provided the regulator 
transforms as a tensor under gauge transformations, 
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[Rk,i 



(13) 



This can be established by using background-covariant 
momenta p — > —iD(A) in the regulators, or p^ — >■ 
r^^-' [0, A] as necessary for optimized flows. For $ = 0, 
Eq. ((ni) entails the flow of Tk [A] , 



dtTk[A] = isTr 
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-dtRk : 



(14) 



which is gauge invariant for regulators obeying Eg . ()13p . 
Note that the flow dH]) is not closed [Mill Hi]: the 
right-hand side depends on F)^ ' '[0, A], whereas the left- 

(2) 

hand side only allows for the computation of F)^ ^ = 



F^ ' ^[0,A]. The necessary approximation for the direct 
use of (fn|) is therefore to set 



F 
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The resulting flow ([14]) is indeed closed and can be solved 
within powerful heat-kernel techniques [3J-[3^- For in- 
stance, this approach predicts the existence of an infrared 
fixed point of the coupling at zero and finite temperature 
[36| . as it is similarly found in Fandau gauge QCD. 

However, within the present context of the order- 
parameter potential, it is crucial to go beyond the 
approximation (jisp used in [37| . as the confinement- 
deconfinement phase transition is rather sensitive to 
the correct mid-momentum and infrared behavior of 
the fluctuation-field propagator F^ ' [0, A]. Therefore, 
a proper distinction between the fluctuation-field and 
background-field dependence of the action is mandatory. 



This gauge has several benefits: first, it projects on (co- 
variantly) transversal degrees of freedom, and second, 
the longitudinal components of Green functions decouple 
from the dynamics of the transversal ones, and thirdly 
it is a fixed point of the flow [33]. As the longitudi- 
nal components are subject to modified Slavnov- Taylor 
identities, this minimizes the truncation error. This is 
also related to a second issue, namely the gauge depen- 
dence of the background-field effective action F^ [A\ . For 
Fandau-DeWitt gauge, this action is identical to the ge- 
ometrical effective action, see e.g. [23, |40, |4l|, which 
is gauge invariant also with respect to quantum gauge 
transformations. There, the fluctuation field agrees with 
a — A — A only in leading order. The background-field 
approach in Fandau-DeWitt gauge can indeed be under- 
stood as the leading order of a manifestly gauge-invariant 
approach to functional RG flows [2^ |41[ . 

We proceed with constructing the key input 
Fi,''°'[0,A], 



First we remark that Fj. ' [0,0](y» ) is sim- 
ply the propagator in Fandau gauge which has been com- 
puted on the lattice [5[ as well as by functional methods 
[i-I|. The fuU RG trajectory F^^'"^[0,0](p2) has been 
computed in [^, [g^ for the optimized regulator -Ropt- 
Now the auxiliary background gauge symmetry comes 
to our aid. It constrains the extension of the Fandau- 



gauge two-point function F), ' [0, 0] to F 
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0, A] as the 



latter has to transform as a tensor under gauge transfor- 
mations. We conclude that 
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where Fp^ denotes the held strength tensor in the ad- 
joined algebra, and the function f{x) is non-singular at 
X — 0. Note that covariantly longitudinal correction 
terms in Eq. (J16p are irrelevant as we are in the Fandau- 
DeWitt gauge. The / terms cannot be obtained from 
the Fandau gauge propagator. They are indeed related 
to higher Fandau-gauge Green functions. 

Next we briefly recall the results for the Fandau gauge 
propagators [3, |^, [23, [30] : the ghost and gluon propaga- 
tors can be parameterized as 

rS^[o,o](p') = p'Za{p')Ut{p)i 

+/^nL(rii, (17) 



where 



IV. BACKGROUND-FIELD FLOWS AND 
LANDAU-GAUGE YANG-MILLS THEORY 
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The crucial ingredient for our studies of the con- 
finement phase transition is the two-point function 
T].' '[0,A]. As it is not the output of the pure 
background-field flow ^^, we have to compute it sepa- 
rately. From now on, we restrict ourselves to the Fandau- 
DeWitt gauge ([TU)) with the gauge parameter set to ^ = 0. 



for the gluon and 



r'l'^\oMp')=p'Zcip')i 



(19) 



for the ghost. For the longitudinal dressing function, we 
have Zl = 1 + C(C)- Hence, it drops out of all diagrams 



beyond one loop. In the deep infrared, the dressing func- 
tions Z^c exhibit a leading momentum behavior 
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(20) 



Landau-gauge Yang-Mills theory admits a one-parameter 
family of infrared solutions consistent with RG invari- 
ance Q, the underlying structure still being subject to 
current research. This family of solutions can be param- 
eterized by an infrared boundary condition for the ghost 
propagator, specifying a value for Zc{p^ = 0). This fact 
is reflected in recent lattice solutions on relatively small 
lattices, [42|, and in the strong-coupling limit [43[, for an 
alternative point of view see [441 . 

For Zc{p'^ ^- 0) — >■ 0, there is a unique scaling solution, 
[45l |46| . The two exponents ka and kc ai'c then related 
by the sum rule arising from a non-renormalization the- 
orem for the ghost-gluon vertex |47| . 



0^ KA + 2KC + 



(21) 



in d dimensional spacetime [3, |4^, \^ . Admissible solu- 
tions are bound to lie in the range kc G [1/2 ,1]. In pure 
non-abelian gauge theories, nc has been computed by a 
variety of methods [^-Q . The precise value depends on 
the IR behavior of the ghost-gluon vertex [4^ . In most 
DSE and FRG computations we are led to {d = 4) 

Kc = 0.59535... and ka = -2kc = -1.1907..., (22) 

being the value for the optimized regulator |8|. The reg- 
ulator dependence in FRG computations leads to a kc 
ran ge o f kq S [0.539, 0.595], see Q; for a specific flow, 
see [30]. These results entail the Kugo-Ojima/Gribov- 
Zwanziger confinement scenario. The gluon is infrared 
screened, the propagator even tends to zero, see (|20 |) . ([22]) . 
whereas the ghost is infrared enhanced. Due to the 
non-renormalization property of the ghost-gluon vertex, 
[4 71 148| , a running coupling can be defined in terms of 
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which runs towards an IR fixed point, see Eq. (j20p . In 
Fig. [1] we show the momentum dependence of the ghost- 
and gluon propagator as obtained from a functional RG 
study in comparison to lattice results p . 

A different type of decoupling solution is found for non- 
vanishing Zc{0): here, the gluon propagator tends to a 
constant in the infrared, p^Za{p^) — >■ m^, for related 
work see e.g. [3,|43]. It should be stressed that the gluon 
propagator then does not correspond to the propagator 
of a massive physical particle, but clear indications for 
positivity violation related to gluon confinement are ob- 
served, j3,[50] ■ Still, the gluon decouples from the dynam- 
ics as does a massive particle. The qualitative infrared 
behavior is then characterized by the exponents 



Though the infrared exponents for the scaling (p2|) and 
decoupling solutions ([M|) deviate from another, the prop- 
agators do only differ in the deep infrared. It has been 
suggested in 19| that the infrared boundary condition is 
directly related to the global part of the gauge fixing, 
hence reflecting different resolutions of the Gribov prob- 
lem. Indeed, the infrared boundary condition has been 
used as a global completion of the gauge fixing in [42| . 

For most parts of the present work, the difference be- 
tween the scaling and the decoupling solution is of minor 
importance. For concrete numerical computations, we 
combine information about the propagators as obtained 
from the lattice as well as functional methods, complet- 
ing these propagators in the deep IR with the scaling 
solution. The latter is actually singled out by the re- 
quirement of global BRST for Landau gauge Yang-Mills 
with standard local BRST invariance. 

The discussion so far applies to Yang-Mills theory at 
zero temperature. Several modifications arise in the pres- 
ence of a thermal bath, such as the immediate replace- 
ment of continuous loop energies po by Matsubara fre- 
quencies po ^ LUn = ^imT in the imaginary-time formal- 
ism. Moreover, the gluon propagator acquires an addi- 
tional component, as the contributions longitudinal and 
transversal to the heat bath become independent. 

In this work, we neglect the finite-temperature mod- 
ifications of the propagators, but work with zero- 
temperature propagators evaluated at the Matsubara fre- 
quencies. In scalar theories it has been shown that 
this approximation already provides a quantitative in- 
sight into the finite-temperature phase structure [5l|. 
First results for finite-temperature gluon- and ghost- 
propagators indeed indicate that the propagators are lit- 
tle modified [53,15^ for Matsubara frequencies 2'KTn with 
\n\ > 2,3. Significant changes have been found for the 
gluon propagator longitudinal to the heat bath which is 
increased compared to the transversal counterpart. We 
stress that the inclusion of the full temperature depen- 
dence of the propagators as well as the order-parameter 
fiuctuations is inevitable for an accurate determination 
of, e.g., the critical exponents or the thermodynamic 
properties of the theory (see, e.g., j54|). 

The presence of finite temperature also takes influence 
on the form of the propagator at finite background field, 
as another field invariant, the Polyakov loop L, exists. 
This adds further terms to the right-hand side of Eq. (|16p . 
Speciafizing to constant background fields A = (Ag), we 
find schematically 



(rf")[o,(Ao)])f, 
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- 1 , and Kc = . 



(24) 



(rr^[0,0](-Z?^));it+i terms, 

(25) 

since the / term in ([TS]) vanishes, F((Ao)) = 0, for 
(Ao) =const. In this work we drop the L contribution in 
(PSj) which is related to the second derivative of the order- 
parameter potential via Nielsen identities [2^, |33- We 
expect that this term affects our results only on a quan- 
titative level, e. g., the quality of critical exponents. For 
example in SU(2) gauge theory, we expect a second-order 
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FIG. 1; Momentum dependence of the gluon (left panel) and ghost (right panel) 2-point functions at vanishing temperature. 
We show the FRG results from Ref. [9| (black solid line) and from lattice simulations from Ref. [5| (red points). 



phase transition. Here, the critical dynamics encoded in 
the critical exponents is sensitive to order-parameter fluc- 
tuations as is well knovifn from studies of scalar 0{N) 
theories. The L terms take a direct influence on the 
spectrum of order-parameter fluctuations, so that we ex- 
pect these terms to be relevant at criticality. Indeed, the 
role of order-parameter fluctuations has been studied in 
Ref. 17| for SU{2) where the correct Z2 critical expo- 
nents have been found. However, the phase transition 
in SU{N) gauge theories {N > 3) is of first order and 
therefore less affected by our approximation of dropping 
the L terms. 

Let us finally stress that our approximations at finite 
temperature do not take any infiuence on our conclu- 
sions about confinement in the zero-temperature limit, 
discussed below. In particular, the background covari- 
antization of the transverse propagator in Eq. (|16|) be- 
comes exact in this limit, representing a first important 
result of our present work. This paves the way for a 
fully consistent low-energy RG analysis of QCD in the 
background-field formalism, see Refs. [5a, [5g| for a study 
of 1-fiavor QCD. 

From Eq. (|16p and the above results for the Landau 
gauge propagators we already conclude that the trun- 
cation (jisp is not working well in the (deep) infrared 
(p -C Aqcd). There we expect a fixed point for the cou- 
pling, (|23p , which entails a constant dressing for the prop- 
agator of the background field, F'^'^': background gauge 
invariance leads to RG invariance of gA, and hence to 



that important aspects of the infrared physics can easily 
be missed by truncations based on ([T5|) . 



z;'^ 



(26) 



which results in a constant dressing. Therefore the back- 
ground field propagator l/Fg ' diverges in the infrared 
whereas the propagator of the dynamical fluctuation field 

a, I/Fq ' is suppressed in the infrared for both scaling 
and decoupling solution. Moreover, the ghost propagator 
is infrared enhanced for the scaling solution in contradis- 
tinction to the one-loop truncation used so far in most 
background field fiows, that is Zq ~ 1 which would only 
be compatible with the decoupling solution. We conclude 



V. CONFINEMENT CRITERION 

Our study of the deconfinement phase transition is 
based on an order parameter related to the Polyakov 
loop variable, see Eq. ([T]). The negative logarithm of the 
Polyakov loop expectation value (L) can be interpreted 
as the free energy of a single static color source in the fun- 
damental representation of the gauge group [57| . In this 
sense, an infinite free energy associated with confinement 
is indicated as (L) — )■ 0, whereas (L) ^ signals decon- 
finement. For gauge groups with a nontrivial center, (L) 
measures whether center symmetry is realized by the en- 
semble under consideration jSTj . As (i) transforms non- 
trivially under center transformations a center-symmetric 
(disordered) ground state automatically ensures (L) = 0, 
whereas deconfinement (i) 7^ is related to the breaking 
of this symmetry, pointing to an ordered phase. 

The background-field formalism used in this work al- 
lows us to fix the fluctuations with respect to Landau- 
DeWitt gauge and simultaneously maintain gauge in- 
variance for the background field A which we relate to 
the Polyakov loop in the following manner: we use the 
Polyakov gauge [22] by gauge-rotating the background 
field into the Cartan subalgebra and imposing SqAq = 0. 
From the knowledge of Aq , the value of the corresponding 
Polyakov loop i[Ao] can immediately be inferred. Once 
the effective action Y\A = A\ is constructed within the 
background-field formalism, the minimum of the action 
represents the expectation value of the fluctuating quan- 
tum field, A = A = (A). The corresponding Polyakov- 
loop value then is L[Aq\ — i[(Ao)]. 

In addition to (L), also L[(Ao)] serves as an order pa- 
rameter for confinement: first, L[{Aq)] is an upper bound 
for the Polyakov loop expectation value due to the Jensen 
inequality, L[(Ao)] > {L[Aq]), and therefore is nonzero in 
the deconfined phase. Second, it has been shown in [l7[ 



that L[(Ao)] vanishes identically in the center-symmetric 
phase where also (L[Ao]) = 0. We conclude that together 
with L[{Aq)] also Aq ~ Aq ~ (Ao) is an order parameter 
for center symmetry and confinement in the Polyakov 
gauge. In the following, we indeed concentrate on the 
effective potential V^((Ao)) for this order parameter. 

Let us now recapitulate the confinement criterion put 
forward in Ref. [ij. This criterion relates the IR behav- 
ior of gluon and ghost 2-point functions to the effective 
potential for the order parameter (Aq), starting from the 
full flow as displayed in Eq. ^. 

The following simplified analytical discussion is based 
on the assumption that the second term in Eq. ([4]) pro- 

portional to ~ StF^ ' (cf. Eq. dS])) is sublcading. This 
term with the explicit k integral resembles the full inte- 

(2) 

grated flow except for the substitution dtRk — > ^tT^ . 

(2) 

In the UV, its subleading role is obvious, since dtT^, is 
of order as whereas dtRk is of order one. In the deep 
infrared such an ordering cannot be found. Nevertheless, 
our full numerical study shows that the term depending 

(2) 

on dtTj, is subleading for a study of the Polyakov loop 
on all scales studied in this work. 

(2) 

Anticipating the subdominance of the StF^ term, we 
study the influence of the first term of Eq. (j3]) on the 
Polyakov-loop potential in the UV and IR regime. (The 
remaining terms are irrelevant for this discussion, as the 
potential is finite and does not require counterterms.) In 
the UV regime {p^ ^ T^), perturbation theory holds and 
the inverse propagators of the longitudinal and transver- 
sal gluons and the ghosts are given by F^ j^^^-^ {p^) = p^- 
In the presence of a constant background field (Aq), 
the momentum is replaced by the background covariant 
derivative, i. e. po — > —iDq. With the parameterization 

/3g(AS) ^2nJ2 T'^r = 2^ E T^^v^l^l v' = 1, (27) 

T"6Cartan T"6Cartan 

the spectrum of the background covariant Laplacian be- 
comes 

p2 ^ spec{-Z?2[(Ao)]} = p2 + (2^r)2(n - |0|zy,)2 , (28) 

where n g Z, and (T"-)'"^ = — i/'^''^ denotes the gener- 
ators of the adjoint representation of the gauge group 
under consideration, vi denotes the eigenvalues of the 
hermitian color matrix occurring in Eq. (|27p . 



i^i = spec{(r''z;'')'"=|i;^ = 1} 



(29) 



and therefore depends on the direction of the unit vec- 
tor u°. The index £ labels these eigenvalues, the number 
of which is equal to the dimension dadj of the adjoint 
representation of the gauge group, £ = 1, . . . , dadj, e.g., 
f^adj = N'^ — 1 for SU{N). For each non- vanishing eigen- 
value vg there exists an eigenvalue —vg. For SU{2), we 
have vg = ±1,0. Equation p8)) reveals that 4>^ = \(l)\v'^ 
denotes a set of compact variables, as an arbitrary shift of 
(j)"' can be mapped back onto a compact domain for 0° by 
a corresponding shift of n, i.e., the Matsubara frequency. 



With these prerequisites, the perturbative limit of the 
effective order-parameter potential y in d > 2 dimensions 
is given by 



V 



UV 



{(f,) (2-d)F(|) i^^cos27rn|(/)|i/f 



J^d 



-d/2 



1=1 n=l 



(30) 



Here, we have dropped a temperature- and field- 
independent constant. The dimensionality of the po- 
tential is determined by the dimension of the Cartan 
(sub) algebra. This perturbative V^^ corresponds to the 
well-known Weiss potential [531 , generalized to d dimen- 
sions [59[. It exhibits maxima at the center-symmetric 
points where L[{Ao)] = (and thus also {L) = 0), imply- 
ing that the perturbative ground state is not confining, 
i.e. (L) 7^ 0. Since the eigenvalues Vi are pairwise identi- 
cal with respect to their absolute values, the Weiss poten- 
tial for a given gauge group can be considered as a super- 
position of SU{2) potentials with different periodicities 
determined by the eigenvalues vi . The eigenvalues can be 
viewed as Fourier frequencies of the order-parameter po- 
tential. We stress that this also holds in non-perturbative 
studies of the Weiss potential. Hence we have 



Vi^t^) 



9$I^su(2)(t^;|0l)- 



(31) 



Next we perform the same analysis in the IR. With 
the paramcterizations ((T7)) and P^ . the dressing func- 
tions Za{p^), Zc'{p'^) are characterized by the power-law 
behavior (PD|) in the deep IR, p^ <C Aq(-,q. Quantita- 
tively, the effective potential V^cj)) is dominantly induced 
by fluctuations with momenta near the temperature scale 
p2 ro (27rr)2. At low temperatures (27rT) < Aqcd, the 
flrst term of Eq. ^ thus induces an effective potential 
which arises dominantly from fluctuations in the deep IR, 
characterized by the exponents ka,c- By coupling the 
fluctuations to the background field, p^ — )- — I3^[(Ao)], 
we obtain the following low-temperature effective poten- 
tial from the power-law behavior of the two-point Green 
functions in the deep IR: 



y^«(0) = 1 + 



(d- 1)ka - 2kc 
d-2 



V^^{<t>). 



Compared to the perturbative Weiss potential ([30)) we 
observe that the effective potential is reversed if 



2kc - (d- 1)ka > d-2. 



(32) 



In this case, the confining center-symmetric points of the 
Weiss potential turn from maxima to minima: the order 
parameter acquires a center-symmetric value, such that 
i[(Ao)] = (L) = 0. We conclude that Eq. (|32]) serves 
as a criterion for quark confinement. Provided that the 

(2) 

term ^ F^ that we dropped for this discussion does 
not modify this result, this criterion is sufficient for the 
occurrence of a center-symmetric confining phase at low 
temperatures. 



Let us discuss this criterion in the hght of the IR solu- 
tions for the propagators available in the Landau gauge. 
For the scaling solution, the IR exponents are related by 
the sum rule (|2ip . simplifying the confinement criterion 
to 



Kc > 



d~3 



(33) 



It is instructive to compare this simple criterion for quark 
confinement with related criteria in the d ~ 4 case: 
Kd=i > 1/4. This criterion includes the Kugo-Ojima 
criterion for color confinement k > as well as the 
Zwanziger horizon condition for the ghost k > 0, which 
are both necessary but not sufficient criteria. The cor- 
responding horizon condition for the gluon, k > 1/2, is 
stronger, since the latter is a sufficient but not a necessary 
condition for the transversal gluons to exhibit positivity 
violation. 

For the decoupling solution with kc — and ka = 
— 1 in (i = 4, the criterion ((32|) is satisfied as well and 
hence, the whole one-parameter family of Landau-gauge 
IR solutions is confining. 

To summarize: in color-confined gauge theories, the 
suppressed gluon and the enhanced (or constant) ghost 
fiuctuations in the IR induce an effective potential for the 
Polyakov loop which corresponds to a center-disordered 
confining ground state; this implies an infinite free energy 
for a single quark and thus relates color confinement to 
quark confinement. 



VI. NUMERICAL RESULTS 

In the following, we present our results for the or- 
der parameter L[{Aq)] as a function of temperature for 
SU{N), Sp{2) and ^(7) Yang-Mills theory. Recall that 
-^[(Ao)] > (L[Ao]) in the deconfined phase. Our computa- 
tion of the effective potential V{{An)) for SU{N), Sp{2) 
and E{7) Yang-Mills theory involves several approxima- 
tions which can easily be improved on, once more pre- 
cise propagator data from the lattice or from functional 
methods is available: first, we employ the same solution 
for the ghost and gluon propagators as obtained from a 
functional RG study [3, [231 for all gauge groups, see also 
Fig. [T] In a first approximation, this can be justified, 
since the propagators are identical in leading order in a 
1/A'^ expansion where N is the number of colors. How- 
ever, even for a small number of colors, it has indeed been 
found on the lattice that SU{2) and SU{3) propagators 
agree within errors [601 • 

As a second approximation, wc do not take a possible 
modification of the functional form of zero-temperature 
and finite-temperature propagators into account. We 
also neglect that the transversal gluon propagator splits 
into independent components longitudinal and transver- 
sal to the heat bath at finite T. Our approximation to 
the finite-temperature propagators corresponds to insert- 
ing Matsubara frequencies into the momentum argument 



of the zero-temperature propagator functions. From an 
RG point of view, this represents the zeroth-order ap- 
proximation to the full temperature-dependent propaga- 
tors. Nevertheless, we expect that this already provides 
a quantitative insight into the finite-temperature phase 
structure for the following reasons: finite-temperature 
modifications of the propagators are expected to occur 
for momentum scales below the temperature scale and, 
more prominently, for temperatures below Tc, as is con- 
firmed by corresponding functional and lattice studies 
[52 |. By contrast, the effective order-parameter poten- 
tial is dominantly built up from momentum modes near 
the scale 2ttT. Therefore, detecting Tc from above, the 
IR properties of the propagators are hardly probed and 
only the decisive mid-momentum region together with 
the perturbative high-momentum tail of the propagators 
effectively enters in the computation of the potential. Ne- 
glecting a potentially strong explicit temperature depen- 
dence of the propagators for small temperatures T <^Tc 
and/or momentum modes below the temperature scale is 
thus an acceptable approximation for detecting the phase 
boundary. The validity of this approximation has been 
verified explicitly for scalar theories in Ref. [51| . 

In addition to our analytical discussion of the con- 
finement criterion, we have now used the full functional 
flow equation including the term depending on dtTl. in 
Eq. ^ in our numerical study. We observe that the order 
of the phase transition for a given gauge theory remains 
unchanged upon the inclusion of this term. Moreover 
the phase transition temperature increases only by < 7% 
when this term is added. For a qualitative understanding 
of the order-parameter potential as discussed in the pre- 
ceding section, the omission of this term is hence justified 
which confirms the picture arising from the our confine- 
ment criterion. For details on the numerical computation 
of the order-parameter potential, wc refer to App. |^ 

In order to convert our results into physical units, we 
fix our propagators relative to the lattice scales. In turn, 
the propagators on the lattice can be converted into phys- 
ical units by measuring lattice momenta in units of, e.g., 
the string tension. In this manner, we can determine 
Tc in physical units corresponding to a string tension of 
a = 440 MeV. In our studies we keep the position of the 
peak of the gluon propagator fixed for all gauge groups. 
This provides a prescription for a comparison of lattice 
results for Tc and our results. 



A. SU{N) 

Let us first consider the gauge groups SU{N). In 
Fig. [21 we show our results for the order parameter 
L[(Ao)] for N = 3,5,7,9,11 as a function of T/Tc, the 
corresponding results for SU{2) and SU{12) can be found 
in Figs. [3] and IH respectively. We find a second order 
phase transition for SU{2) and a first-order phase transi- 
tion for SU{N) (TV = 3, 4, . . . , 12). For SU{2) the phase 
transition occurs at Tc ~ 265 MeV. For S'C/(3) we find 
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FIG. 2: Polyakov loop L[(Ao)] as a function of temperature 
for SU(3), SU(5), SU(7), SU(9), SU(ll). We observe that the 
phase transition is of first order. Hardly any difference for 
the order parameter is visible for N > 5, suggesting a close 
proximity of these gauge groups to the large- A'' limit. 
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FIG. 3: Polyakov loop L[{Ao}] for SU{2) (blue/dashed line) 
and Sp{2) (black/solid line). The phase transition is of second 
order for SU{2) and of first order for Sp{2). 



Tc ~ 291 MeV. With increasing rank of the gauge group, 
the phase transition temperature increases slightly and 
approaches T^ sa 295 MeV for SU{h). For iV > 5 our 
results for the order-parameter are essentially indepen- 
dent of N. In other words, our results for the phase 
transition temperature for SU{2>) is already close to the 
large- A^ value. This independence of the Polyakov loop 
on A^ for A^ > 5 is in accordance with recent lattice stud- 
ies of SU{N) Yang-Mills theories [61[. Since we employ 
the same propagators for all SU{N), the increase of Tc is 
only due to the increase in the rank of the gauge group. 
In accordance with the weak dependence of Tc on A^ > 2, 
the order-parameter as a function of T/Tc depends only 
slightly on the rank of the gauge group. Note that our 
result for the Polyakov loop T[(Ao)] for T/Tc > 1 is 
higher than the corresponding expectation value (T[Ao]) 
of the Polyakov loop as obtained from lattice simulations, 
being in perfect agreement with the Jensen inequality 
L[{Ao)] > (T[Ao]). 

At this point we would like to emphasize that our stud- 
ies are of course not bound to A^ < 12. Our approach can 
be straightforwardly generalized to iV > 12 with the aid 
of Eq. (|3ip . Our numerical study of a given gauge group 
involves three simple steps: computing Vgu(2) : finding the 
eigenvalues of the generators of the Cartan subalgebra in 
the adjoint representation for the gauge group under con- 
sideration, and finally minimizing Eq. (j3ip . Therefore the 
computation of the order parameter for very large gauge 
groups is not considerably more involved than for smaller 
ones. 

In view of the approximations listed above, we ex- 
pect corrections to our results from modifications of the 
propagators due to finite temperature and due to order- 
parameter fluctuations. Whereas finite-temperature cor- 
rections of the propagators affect the results for the order 
parameter of all gauge groups, order-parameter fluctua- 



tions play a particularly important role in SU{2) since it 
has a second-order phase transition. 



B. Sp{2) 

The fact that confinement and center symmetry are re- 
lated naively suggests that gauge groups with the same 
center may show similar phase transition properties. This 
is, however, not the case as the prime counter-example 
of SU{2) vs. the symplectic group 5*^(2) demonstrates: 
both gauge groups have the same center Z(2), but ex- 
hibit qualitatively different phase-transition properties^. 
Our results for the order parameter L[(Ao)] as a func- 
tion of T/Tc for SU{2) and Sp{2) are depicted in Fig. [3 
The generators of Sp{2) are given in App. |B] We find a 
second-order phase transition for SU{2) and a first-order 
phase transition for Sp{2). Therefore, SU{2) falls into 
the Ising universality class [53| but Sp{2) does not. More- 
over the phase transition temperature for Sp{2) gauge 
theory is close to the value of 5C/(3) gauge theory; we 
find Tc « 286 MeV. Since we use the same propaga- 
tors for both gauge groups and the center of both groups 
is Z{2), it is natural to relate the different nature of the 
phase transition to the different dimensionality of the two 
groups [211 . In fact, the number of degrees of freedom 
in the deconfined phase is much larger in Sp{2) than in 
SU{2). This strong mismatch in the number of dynam- 
ical degrees of freedom in the confined and deconfined 
phase appears to enforce a first-order phase transition in 
Sp{2). In this respect these findings resemble the situa- 
tion in the case of SU{3) in 2 -f- 1 and 3-1-1 space-time 
dimensions. While the phase transition in SU{3) is of 



^ In our conventions Sp{l) is isomorphic to SU{2). 
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FIG. 4: Polyakov loop L[(Ao)] for S'[7(12) (blue/dashed line) 
and E{7) (black/solid line). The phase transition is of first 
order for both SU{12) and E{7) gauge theory. 



first order in 3 + 1 dimensions, it is of second order in 
2 + 1 dimensions [62[- Again, this might be traced back 
to the fact that the mismatch in the number of dynami- 
cal degrees of freedom in the deconfined phase is smaller 
in d = 2 + 1 than it is in d = 3 + 1. 

Our findings for the nature of the Sp(2) phase transi- 
tion are in accordance with lattice simulations [20|, UM, ■ 
Quantitatively, our approximation of neglecting terms 
oc y'((Ao)) on the right-hand of Eq. Q, which account 
for order-parameter fluctuations, might be more severe 
in Sp{2) due to its similarity to SU{2). The inclusion of 
these fluctuations may only lead to a weaker first-order 
jump and an increase of the critical temperature. 



C. E{7) 

Another interesting test of the proposal that the order 
of the phase transition is related to the size of the gauge 
group [2l| is the following comparison between SU{12) 
and E(7) gauge theory. The dimension of the adjoint 
representation of these two gauge groups is about the 
same: We have dadj = 133 for E{7) and dadj = 143 for 
SU (12). The gauge groups differ with respect to their 
center, being Z{2) for E{7) and Z(12) for SU{12). 

In Fig. 131 our result for the Polyakov loop L[{Ao)] for 
both theories is depicted as a function of T/Tc- As 
discussed above, the phase transition in SU{12) is of 
first order and occurs at Tc « 295 MeV. For E{7), our 
RG approach predicts a first-order phase transition at 
Tc « 295 MeV as weU. 

Our study is thus compatible with the suggestive rela- 
tion of the order of the phase transition and the mismatch 
of the number of degrees of freedom above and below the 
phase transition - provided the glueball spectrum below 
the phase transition in E{7) is similar to that of SU{N). 

However, we also observe that the height of the jump 



of the order parameter is smaller in E{7) than in SU{N) 
for all values of N studied in the present paper. Even 
though the Polyakov loop I/[(Ao)] is not an RG invariant 
quantity, our approach of studying the associated eigen- 
value distribution allows us to give the height of the jump 
a physical meaning. This suggests that the mismatch in 
the number of degrees of freedom is not the only mecha- 
nism that determines the nature of the phase transition. 

In order to gain a better understanding of the nature 
of the phase transition, we study the eigenvalue distribu- 
tion iV(|i/^|) of the spectrum of the color matrix in the 
Cartan subgroup as defined in Eq. ([29]) . In Fig. [5] we 
show A^d^-fl) as a function of the normalized eigenvalues 
Wt/i^imaxl at the ground state of the order-parameter po- 
tential for T — > T^ , approaching the critical temperature 
from above. Here, the eigenvalues have been binned with 
a bin size of Ai'f = 0.005. 

The eigenvalues correspond to Fourier frequencies of 
SU(2) Weiss potentials, cf. Eq. jSOl) and Eq. §^. For 
instance, the higher the dominating Fourier frequency in 
the high-temperature phase, the closer the minimum is 
to (/) = 0, implying that L[(Ao)] is closer to L[{Aq)] = 1. 
By contrast, if lower eigenvalues dominate, L[(Ao)] can 
approach the center ordered state in a smoother fash- 
ion. This is precisely what we observe for E(7) in con- 
tradistinction to the SU{N > 2) gauge groups, where the 
eigenvalues cluster around fe/fimax ~ 0,0.25,0.5,0.75, 
leading to an almost constructive interference of SU{2) 
potentials with almost identical periodicity. We stress 
that the eigenvalue distribution A^di/fj) depends on the 
actual position (Ao)min of the ground-state of the poten- 
tial. Therefore the eigenvalue distribution and hence the 
strength of the first-order phase transition depends on 
the actual trajectory (Ao)min(r) of the physical ground- 
state close to Tc in the space spanned by the generators 
of the Cartan subalgebra. Of course, since L[{Ao)] pro- 
vides only an upper bound for (L) it is not immediately 
clear whether the difference in L[(Ao)] which we observe 
for E{7) and SU{N) also translates into a similar dif- 
ference in (L). If so, we expect the phase transition for 
E{7) to be smoother than for SU{N). Taking into ac- 
count that order parameter fluctuations dropped so far 
can smoothen the phase transition even further, our re- 
sults may not even be taken as a strict excluding evidence 
for a second order phase transition in E(7). 

With respect to our study of Sp{2) we indeed find 
that the eigenvalue distribution 7V(|i/£|) exhibits a pat- 
tern very similar to the one of SU{3) resulting in a jump 
of the order parameter at the phase transition with a 
height comparable to the one of SU{3) Yang-Mills the- 
ory, see Figs. [2] and [3l 



VII. CONCLUSIONS 

In the present paper we have discussed the nature of 
the phase transition in various gauge groups and based 
on a simple confinement criterion put forward in Ref. [l| . 
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FIG. 5: Eigenvalue distribution Af(|i/f|) of the spectrum (pOJ) 
as a function of the (normalized) eigenvalue Ive/vimaxl for 
E{7) (lower panel) and SU{12) (upper panel) at the ground 
state of the potential for T — >■ T^ . 



For our study of the critical dynamics of non-abelian 
gauge theories we computed the order-parameter poten- 
tial in Landau-DeWitt gauge by employing gauge corre- 
lation functions. In particular, for the question whether 
the system is in the confining phase near zero temper- 
ature, the quantitative knowledge of zero-temperature 
propagators is sufficient. Even for questions related to 
the nature and quantitative properties of the phase tran- 
sition near the critical temperature, we have argued that 
we expect that zero-temperature propagators provide for 
a reasonable approximation. Of course, the inclusion of 
knowledge about thermal propagators in our formalism 
is straightforwardly possible. Moreover we have dropped 
order-parameter fiuctuations which, though irrelevant for 
the question of confinement near zero temperature, we 
expect to affect, for instance the quality of, e. g., critical 
exponents. The fiuctuations have been included in the 
S'C/(2)-study in Polyakov gauge, [131, ^'^'^ I'^^d ^^ Ising- 
class critical exponents as expected. 

In agreement with lattice simulations we have found 
a first-order phase transition for SU{N) gauge theories 
with N > 3. For SU{2) the phase transition is of second 
order. Moreover we have studied 5*^(2) and E{7) gauge 
theory and compared our results to SU{N) Yaiig- Mills 
theory. In agreement with lattice simulations [20, l21| 
we observe a first-order phase transition in 5*^(2). As a 
new set of characteristic quantities of the phase transi- 
tion, we have introduced the distribution of eigenvalues 
of (Aq) which within the Polyakov gauge can be related to 
an order parameter. These eigenvalues serve as Fourier 
frequencies of a superposition of SU{2) Weiss potentials 
yielding the full nonperturbative Weiss potential V^[(Ao)] 
from which L[(Ao)] can be deduced as an upper bound to 
the Polyakov loop order parameter. Whereas the order- 
disorder nature of the dcconfinement transition is related 
to center symmetry, the center degrees of freedom them- 
selves are not always the relevant degrees of freedom to 



understand the phase transition. For instance, this is ob- 
vious in the case of Sp{2) which has a center Z(2) but 
exhibits a phase transition different from SU{2) with the 
same center. This is in fact illustrated by the eigenvalue 
distribution: for Sp{2) the eigenvalue distribution at the 
phase transition is similar to the one of SU{3) as are 
other properties of the phase transition. In this picture 
the order-parameter potential can be considered as a de- 
structive interference/superposition of SU{2) potentials 
favoring a first-order phase transition. Moreover, we have 
a stronger mismatch in the number of the dynamical de- 
grees of freedom in Sp{2) in the confined and deconfined 
phase compared to SU{2) Yang-Mills theory f20l.[2H. 

For E{7) gauge theory we find that the phase tran- 
sition is of first order as well. Here, the mismatch in 
the number of dynamical degrees of freedom in the con- 
fined and deconfined phase is even stronger than it is 
in Sp{2), suggesting that the first-order phase transition 
is even stronger. However, our RG study suggests that 
the first-order phase transition is weaker for £'(7) than it 
is for SU [12) or 5*^(2). We have argued that this weak 
first-order transition can be traced back to the eigenvalue 
distribution at the phase transition. In contrast to Sp(2) 
and SU{N) we have found that the distribution exhibits 
distinct equidistant maxima resulting in an almost con- 
structive interference of SU{2) potentials. In this respect 
E{7) is closer to SU{2) than to SU{N) with iV > 3. How- 
ever, further studies arc needed to establish this picture. 
For example, a RG study of SU{3) and 5*^(2) Yang-Mills 
theory in 2 -|- 1 dimensions may help to shed more light on 
the underlying mechanisms of the deconfincmcnt phase 
transition since it is known from lattice simulations that 
the nature of the phase transition in both gauge groups 
changes from first to second order when the number of 
dimensions is reduced [2G, EH E3 . 



Another interesting case is the gauge group G(2) with 
non-trivial center. In this case it has been found [20, IS^ 
that the Polyakov loop exhibits a jump but is non- 
vanishing for all temperatures. A verification of our 
quark confinement criterion with the aid of G(2) Yang- 
Mills theory is under way and will help us to establish 
our findings and to improve our understanding of con- 
finement in gauge theories. 
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Appendix A: Details on the computation of the 
order-parameter potential 



In this addendum we discuss some details on the com- 
putation of the order-parameter potential. The order- 
parameter potential can be obtained directly from an 
evaluation of Eq. ^. While the first term in Eq. ^ 
is independent of our choice of the regulator function, 
the second term is not. Since we employed the opti- 
mized regulator (|5]), we encounter expressions involving 
unit-step functions in the second term on the right-hand 
side of Eq. Q . These unit-step functions depend on the 
background field (Aq). This dependence on (Aq) gener- 
ates divergent RG flows for A — ?> cxd. In principle one can 
deal with these divergences by computing appropriate 
counter-terms at the initial UV scale. In the present pa- 
per we sought for a different approach to circumvent this 
problem and introduced 'smeared' unit-step functions: 



fe{x[{Ao)le) = e-(-[<'^«>l)\ 
lim/0(x[(Ao)],e) = 9(1 - x[{Ao)]) , 



0.207 ^ 



0.2065 



^ 0.206 




0.2055 



0.205 



FIG. 6: <^min as function of the 'smearing' parameter e for 
T = 300 MeV. The numerical data is depicted by dots. The 
red (solid) and blue (dashed) line are the results from the fit 
to the functions given in the first and second line of Eq. (|A2|) , 
respectively. 



(Al) 



where a;[(Ao)] is an arbitrary function depending on the 
background field (Aq). Using fe{x[{Ao)],e) instead of 
0(1 — a;[(Ao)]) yields an order-parameter potential pe- 
riodic in (Aq) for any finite value of e and allows to 
get conveniently rid of the unphysical divergent parts of 
the flow. For our numerical study of the deconfinement 
phase transition we have used e = 7. In Fig. [5] we il- 
lustrate the dependence of the position of the minimum 
'/'mill = /?(Ao)min/(27r) of thc potential on the 'smearing' 
parameter e for T = 300 McV for SU{2) Yang-MiUs the- 
ory. From an extrapolation of our results to e ^> oo using 
the two functions 



C, = 



C. = 



/ 10 

-10 



V 



i 
-i 





'^min 



(e) = 






const. + ae 
const. + - 



-be 



(A2) 



where a, b and c are fit parameters, we estimate that the 
theoretical error is less than 1% when e = 7 is used. Note 
that the fit function in thc second line of Eq. (jA2p can 
be formally deduced from a Taylor expansion around e of 
the integral of a general polynomial in a;[(Ao)] weighted 
hy fe{x[{Ao)],e). 



Ci 



C. 







i 

i 





1 
0-10. 



Appendix B: Generators of Sp{2) 

Our definition of the generators of Sp{2) in the funda- 
mental representation is as follows: 



Cfi = 







i 

-i 



Ci = 



i 

i 







Ct 



I ' 
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V 



'V2 



71 \ 










/ 
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